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Abstract 



We re-examine the issue of the pQCD factorization scale in the exclusive p production 
| in diffractive DIS from the /^-factorization point of view. We find that this scale differs 

significantly from, and it possesses much flatter Q 2 behavior than widely used value 
(Q 2 + m 2 p )/A. With these results in mind, we discuss the Q 2 shape of the p meson 
production cross section. We introduce rescaled cross sections, which might provide 
further insight into the dynamics of p production. We also comment on the recent ZEUS 
q . observation of energy-independent ratio 0"(7*p — > PP)/o'tot( / y*p)- 

i : 1 Introduction 

Ph. 
D . 

^ c| The exclusive production of vector mesons in diffractive DIS 



ry(*)p — » Yp 

turned out to be an ideal testing ground pfl |H] of many predictions of the famous 
color-dipole approach, Within this formalism, the basic quantity is the cross 

section of the color dipole interaction with the target proton o"^ p (x,r), which can be 

approximately related to the conventional gluon density G(x g ,Q ) inside the proton, 
[I]. The hard scale Q , at which the gluon density should be taken, is set both by the 
virtuality and the mass of the vector meson, and in the case of heavy quankonium is 
approximately equal to \{Q 2 + m v ). The latter result can be also obtained in the direct 
DGLAP-inspired calculations of the relevant Feynman diagrams, see, for example, the 
calculations of the exclusive production of J/tp mesons in diffractive DIS [7j. 

Although the identification of the pQCD factoriozation scale with \{Q 2 + rriy) is 
valid only in the case of non-relativistic vector mesons, there is a hope that the same 
relation will hold for the light vector mesons and for large virtualities as well. This hope 
leads to the prediction of a remarkable universality: the cross section of light and heavy 
mesons will exhibit the same Q 2 behavior, if plotted against Q 2 + my rather than Q 2 
alone. 
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On the experimental side, the first data showed that cross sections of p, u>, 4> an d 
J/ij) mesons, taken at equal values of Q 2 + m v and corrected by corresponding SU(4) 
factors, were indeed very close to each other, It is interesting to note that a similar 
dependences on Q 2 + m v were observed not only in the magnitudes of the cross sections 
of different vector meson production, but also in the values of the energy growth expo- 
nent, which is related to the Pomeron intercept, as well as in the slopes of the forward 
diffraction cone. 

However, with the advent of new, more accurate data of p and J/ip production, it 
became clear that this scaling was only approximate. Recently, ZEUS concluded [2] that 
the J/ip production cross section is typically 40% higher that the cross sections of the 
light vector meson production. This difference is especially obvious on the plots of &l 
and gt separately, |Hj. 

A natural question arises as to what extent the universality of the Q 2 = (Q 2 + my)/4 
as the relevant hard scale is valid. Since there is almost no room for doubt in the 
case of J /if) production, the question can be formulated as what is the relevant pQCD 
factorization scale of the high Q 2 production of p mesons . 

The answer to this question was, in fact, given already in [3]. Within the color dipole 
formalism, the pQCD factorization scales for the longitudinally and transverse produced 
p mesons were shown to be Q 2 (pl) ~ 0.15- (Q 2 + m v ) and Q 2 (pt) ~ (0.07 -=-0.1) • (Q 2 + 
m v ). This is notably smaller than (Q 2 + my) /A, and suggests that even the highest Q 2 
experimental points in p production correspond to, at most, semiperturbative values of 



In fact, the exact value of the pQCD factorization scale can be affected by the shape 
of the unintegrated gluon distribution. Unfortunately, when 3 ; appeared, no numerically 
reliable parametrizations of dipole cross section or of the unintegrated gluon structure 
function were available, and one was bound to a semiquantitative guess. The situation 
changed two years ago, when numerially accurate, simple, and ready-to-use parametriza- 
tions of the unintegrated gluon structure function T{x g , k) were obtained from the anal- 
ysis of proton structure function i*2 P , ^01- These parametrizations were devised for 
xsj < 0.01 and for the entire domain of relevant Q 2 values. They were put in the basis 
of the /^-factorization calculations of both light and heavy vector meson production cross 
sections and yielded rather good description of the available data, ^2 • These 

fits now allow for a quantitative reanalysis of the hard scale in the p meson production. 
This is performed in the present paper. 

Numerically accurate understanding of the hard scale in the p production is demanded 
in many applications. For example, a need for the hard scale arises when one attempts 
to understand the Q 2 behavior of the p-meson electroproduction. The early data could 
be successfully parametrized by a simple law 



with n = 2.32 ± 0.10 (ZEUS, [Hj) or n = 2.24 ± 0.09 (HI, [T5]V However, with the 
advent of high precision data in a much broader Q 2 region, it became clear that such 
powerlike fits have very limited applicability domain. A natural question has been raised 

x To be definite, we will deal with p mesons, but the general conclusions are obviously valid for all light 
vector mesons. 
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as what would be the most insightful and physically motivated fit to the Q 2 behavior of 
cr(7*p -» pp). 

In virtually any description of the exclusive p meson production, one has to deal with 
the gluon content of the proton. Clearly, in order to experimenally study the underlying 
mechanism of the vector meson production, one might want to get rid of the "trivial" 
source of the Q 2 behavior that arises due to the gluon density. A simple way to do 
this would be to divide the cross sections by the conventional gluon density squared and 
study the Q 2 behavior of the result. However, this procedure requires knowledge of the 
hard scale, at which the gluon density should be calculated. 

Another issue that demands the understanding of the hard scale in p production is 
related to the energy behavior of the ratio 

_ cr( 7 (*)p -> pp) 

The pQCD improved Regge model predicts that these cross sections rise with energy as 

a{l (*) p ^ ^ oc (W 2 ) 2A *> ; <7 tot (7 W p) oc (W 2 ) A v , (2) 

where Aip = ajp — 1 is the Pomeron intercept, so that the ratio is predicted to rise as 

r p oc (W 2 ) A *> . (3) 

In last years, a large amount of data became available on the p meson photo- and 
electroproduction in a broad range of the total photon-proton energy W and the photon 
virtuality Q 2 , [141115] . This enabled the ZEUS collaboration to study the energy behavior 
of r p at several values of Q 2 . The results, posed an apparent challenge: throughout 
the whole Q range, from Q 2 = up to Q 2 = 27 GeV 2 , the ratio r p was found to be 
energy independent. The conclusion of [S] was that neither Regge model, nor pQCD 
approach can explain this constancy. 

However, the Pomeron intercept is known to significantly depend on the hard scale 
involved in the interaction, see experimental data |14[ I15| I16( I17j and results of the 
phenomenological analysis of 10 j. Thus, when testing prediction (j3J), one must make 
sure that both intercepts in @ are taken at the same hard scale. 

In [0] these scales were identified with Q 2 for the total virtual photoabsorption cross 
section and with Q p = Q = (Q 2 + m p )/A for the p meson production cross section, so 
that the experimental results were given for the quantity: 

= <r(7(*)p-> pp)(W 2 ,Q 2 ) 
<?tot{l^ P ){W 2 ,Q 2 ) 

Note that here the p meson production cross section and the total virtual photoabsorption 
cross section are taken at different virtualities. 

As we argue in this paper, the true scale of the p production can noticeably deviate 
from (Q 2 + m 2 )/4:, especially at very small and very large Q 2 . Thus, it appears that the 
mismatch of the scales can be at least one of the sources of the discrepancy observed. 

The structure of the paper is the following. In Section 2 we briefly review the results of 
the fcf-factorization approach to the exclusive production of vector mesons in diffractive 
DIS. For future guidance, we also show how these results simplify in the case of heavy 
vector mesons. In Section 3 we study which gluons are relevant for the interaction, and 
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settle the factorization scales Q L and Q T for longitudinal and transverse amplitudes for 
any Q 2 . Section 4 serves as an application of the results obtained to the study of the Q 2 
behavior of the p meson electroproduction cross section. In this Section we also propose 
new quantites, the rescaled cross sections, which might provide additional insight into 
the mechanism of the interaction. Finally, Section 5 is devoted to the discussion of the 
results and to the conslusions. 

2 Exclusive production of vector mesons 

Within the familiar color dipole formalism 01 E] the production of a vector meson is 
viewed as a three-step process. First, at distances of the order of the coherence length 
£ c ~ l/(rriNXBj) (where xsj is Bjorken x and is the nucleon mass) upstream the 
target, the incident photon splits into a qq pair, a color dipole. Then, interaction between 
the color dipole and the target takes place. This qq pair receives longitudinal momentum 
transfer, so that its invariant mass squared turns positive and close to the mass of a 
corresponding vector meson. Then, at large formation distances, this qq pair projects 
onto hadronic final states. 

At high energies, the typical longitudinaly size of the interaction region, of the order 
of radius of the proton, is much smaller than both the coherence length and the formation 
length. This allows one to treat the color dipole frozen during the interaction, and the 
basic quantity that appears in the approach is the diagonal color dipole cross section 
(Jdip{x^r)^ which is related to the gluon density of the target proton. The analysis of [HI 
yielded the result that the pQCD factorization scales for the longitudinally and transverse 
produced p mesons were Q {Pl) « 0.15-(Q 2 + m^) and Q (p T ) « (0.07-0.1) • {Q 2 + m 2 v ). 

In this paper we reexamine this issue within the fet-factorization approach, which is 
closely related to the color dipole formalism. As said above, our analysis is based on the 
parametrizations of the unintegrated gluon structure function that were given in jlOj . 

2.1 The basic amplitude 




p p p p 



a) b) 

Figure 1: The QCD-inspired diagrams for 7*p — > Vp process with two gluon t-channel. Only 
Diagr.(a) contributes to the imaginary part of amplitudes. 

Within the fcj-factorization approach, the imaginary part of the amplitude of the 
vector meson production is given by discontinuity of the diagram shown in Fig. The 
upper blob should be understood as sum of four diagrams, shown in Fig. [21 We denote 
the quark and gluon loop transverse momenta and the momentum transfer by k, k, and 
A, respectively. The fraction of the photon lightcone momentum carried by quark is 
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Figure 2: The content of the upper blob in FigOJi in the pQCD approach. 



denoted by z, while the fractions of the proton light cone momentum carried by the two 
gluons are x\ and x 2 . With this notation, the imaginary part of the amplitude can be 
written in a compact form: 



ImA = s 



-^-a s (q )J r (xi,x 2 ,ki,K2) 



dzd 2 k 



i>v( z ->k) ■ J(Ay, A 7 ) . (4) 



4vr 2 J k 4 ~ ovl J z (l-z) 

where k\ j2 = k± ^A. The helicity-dependent integrands I(Ay,A 7 ) have form 

il-2z) 2 2m 



I S (L,L) 
I s (T, T) 

I S (L,T) 

I S (T,L) 
where 

$2 
*1 



4QMz 2 (l - z) 



1 + 



4z(l - z) M + 2m 



<h; 



(eV*)[m 2 3> 2 + + (1 - 2z) 2 (kV*)(e$i 



M 



M + 2m 



-(efc)(V*$i) + 



2m 



M + 2m 



(k£)(kV*)$ 2 ; 



2Mz(\ - z)(l - 2z)(e$i) 



1 + 



(1 - 2z) 2 2m 
4z(l - z) M + 2m 



Mm 
M + 2m 



-2Qz{l - z)(l - 2z)(V*k) 



M 



M + 2m 



$2 • 



;i-2«)(e!fe)*2i 



(5) 



+ 



(r + K) 2 + e 2 (r-£) 2 + e 2 (f+A/2) 2 + e 2 (r-A/2) 2 + e 2 ' 



r + k 



r — k 



r + A/2 r-A/2 



'+k) 2 + £ 2 (r-£) 2 + e 2 (f+A/2) 2 + e 2 (r-A/2) 2 + e 2 ' 



with r = k — (l — 2z)A/2 and e 2 = z(l — z)Q 2 + m 2 . Finally, the strong coupling constant 
is taken at q 2 = max[e 2 + k 2 , k 2 ]. 

In the absence of A — k correlations, and for a very asymmetric gluon pair, the off- 
forward gluon structure function T (x%, x 2 , K 2 ) that enters (jU can be approximately 
related to the forward gluon density T{x 9 , it) via 

T{x\,x 2 < xi,ki,k 2 ) ~ F(x gi it) ■ exp (--bsypfc 
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Here x g ~ 0.41xi; the coefficient 0.41 is just a convenient representation of the off- forward 
to forward gluon structure function relation found in |18j . Numerical parametrizations 
of the forward unintegrated gluon density F(x g , k) for any practical values of x g and k 2 
can be found in ^Uj- The slope 631P can be experimentally accessed in the high-mass 
elastic diffraction. 

The vector meson wave function ipv(z, k) describes the projection of the qq pair onto 
the physical vector meson. It is normalized so that the forward value of the vector meson 
formfactor is unity, and the free parameters are chosen to reproduce the experimentally 
observed value of the vector meson electronic decay width r(V — > e + e~). In what 
concerns the shape of the radial wave function, we followed a pragmatic strategy. Namely, 
we took a simple Ansatz for the wave function, namely, the oscillator type wave function 
and performed all calculations with it. In order to control the level of uncertainty, 
introduced by the particular choice of the wave function, we redid calculations with 
another wave function Ansatz, namely, the Coulomb wave function, and compared the 
results. Since these two wave functions represent the two extremes (very compact and 
very broad wave functions that still lead to the same value of the electronic decay width) , 
this difference gives a reliable estimate of the uncertainty. It is given typically by factor 
of 1.5 for magnitudes of the cross sections, while in the observables that involve ratios of 
the cross sections (including slopes and intercepts) this uncertainty is reduced. Details 
can be found in |13[ II 1 j . 

Note also that when deriving ©, we treated the vector mesons as 15 wave states 
and used the pure S- wave spinorial structure instead of 7^, |19j . 



2.2 The heavy meson limit and the scaling property 

It is instructive to look at the production amplitude @ in the heavy quarkonium limit. 



In this case, the vector mesons is essentially non-relativistic, z — > 1/2 , k 2 <C fh 2 and 



e 2 = z(l - z)Q 2 + m 2 « ^(Q 2 + m 2 v ) = Q 2 . 

Even in the photoproduction limit, the mass of the meson sets the hard energy scale of the 

— 2 

interaction, and one should expect the DGLAP-inspired leading log Q approximation 
to be a viable approach to the problem. The DGLAP approximation of the heavy vector 
meson production implies the following ordering 

P«k 2 «Q 2 , 

which leads to simplifications in the production amplitude. Consider the forward scat- 
tering. In this case only helicity conserving amplitudes survive, and one obtains 



?2 8Qm v 2 8m v 

* {Q 2 + m 2 v ) 2 ' ( ' (Q 2 + m 2 ^ 



We reiterate that this simplification works while k 2 is sufficiently small, namely, k 2 <C Q 2 ■ 
In the spirit of the DGLAP approach, the integrands I(Ay,A 7 ) do not depend on k 2 , 
aside from the overall R 2 factor, whose origin is the color neutrality of the proton. Then, 
integration of k 2 immediately leads to the conventional gluon density: 

^a s (K 2 )^(x 9 ,K 2 )/(T,T) -> ( g 2 8 ^^ 2 J q ^a s (K 2 )^(x g ,K 2 ) 
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The integration of the wave function can be related to the electronic decay widths of the 
vector meson. The result for the forward differential cross section reads, [7]: 



da 

W\ 



7T 



t=0 



-my T(V — ► e + e 



a s (Q z )-G(x g ,Q z ) 



Q 6 



(6) 



Note that in this limit the difference between the pure 5-wave vector meson and the 
vector meson with 7^ coupling vanishes. 

Eq.© identifies the hard scale that defined the gluon content of the proton in such 
a reaction with the quantity Q = (Q 2 + m 2 ) /4. Although this derivation is valid only 
for very heavy vector mesons, it is often hoped that for the light mesons the hard scale 
will have the same form. We show below that this is not true. 

This result, although an approximation, suggests the scaling property of the vector 
mesons production: the Q 2 behavior of all vector mesons cross sections should be pro- 
portional to each other if taken at equal Q 2 + m v . Note that © does not suggest that 
the absolute values of the cross sections should be equal, even when corrected by the 



relevant flavor content factors c 
predict universality of my ■ T(V 



V 



simply due to the fact that this approach does not 
e+e'). 



3 The scale in exclusive p production 

In this section we perform an analysis of the amplitude Ijljl. taken for simplicity at A = 0, 
and find what values of the transverse gluon momenta are essential in the interaction. 
These values will define the hard scale of the conventional gluon density. Again, this 
analysis is essentially similar to what was done in [3]. 



3.1 Weight functions: mapping the glue 

We start with a qualitative analysis of the production of longitudinal p mesons at high 
Q 2 . Let us for a moment neglect the scaling violations and the other sources of marginal 
k 2 behavior, and rewrite the amplitude @ in the following simplified form: 

dZ 2 _ .„ 1 r 1 



/dR If /"-»-» 

— F(x g , k) ■ W l (k 2 ) ; W L ^^ dzz(l-z) d 2 k^{z, k) $ 2 
K K JO J 



with 



$2 = ~ : r - — : r + 



(k + k) 2 + e 2 (k- k) 2 + e 2 k 2 + e 2 

At large Q 2 , the most of the k 2 integration comes from the logarithmic region with 
k 2 -C k 2 <C s 2 . This leads to simplifications: 
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W l (k 2 )^ I dzz(\-z) I d 2 W(z,k)- " + (7) 



If we denote by (e 2 ) the typical values of e 2 = z(l — z)Q 2 + m 2 that dominate the 
integral (JJJ), then the weight factor Wl(k 2 ) stays almost constant at k 2 <C (e 2 ), and 
quickly decreases for k 2 > (e 2 ). It should have a form of "smoothed step function", and 
effectively cuts off from above the k 2 region essential for the interaction. If we further 
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approximate it by the exact step function 9(k 2 — Q L ) with properly defined cut-off scale 

2 

then the integration can be done exactly: 

— 2 

/ — Hx g ^)W L (K 2 ) => W L (0) • / _^(x 9 ,K) = iy L (0)-G(x 9 ,Qi). (8) 

JO K Jo Hi 

This transition is precisely what is implied in the DGLAP approach, recall derivation of 
©• Thus, a procedure that would yield the values of the relevant hard scale, at which 
the conventional gluon density should be taken, is the following: take the weight function 
Wl{k 2 ) and find the value Q L that would lead to the most justified replacement ©• 

Obviously, the scale Q 2 L should be close to (e 2 + k 2 ). Due to the broad wave function 
of the p, the typical values of z(l — z) will be less than 1/4. As a result, the scale of the 
gluon density in © should be noticeably softer than (Q 2 + m 2 )/4. 

In the case of the production of the transversely polarized p mesons, the contribution 
of the very asymmetric qq pairs is enhanced. This leads to typical values of z still smaller, 
and the departure of the corresponding factorization scale Q T from (Q 2 + m 2 )/4 becomes 
even more prominent. 



3.2 Numerical results 

The above expectations are corroborated by the numerical analysis. The functions 
Wl{k 2 ) and Wt(k 2 ) are defined now via 

1 I" dH 2 1 f d,K 2 

-ImA L ^ L = / — F(x g , k) ■ W l (k 2 ) ; -ImA T ^ T = / — H%g, ■ W t (k 2 ) 

S J K S J K 

for all values of Q 2 . Their behavior for Q 2 = 100 GeV 2 is illustrated in Fig. |3J where we 
show the ratios Wl(k 2 ) /Wl(0) and Wt{k 2 )/Wt(0) as a functions of k 2 . One sees that 
these weight functions start decreasing at k 2 <C Q 2 and reach 1/2 at k 2 ~ 15 GeV 2 and 
k 2 ~ 10 GeV 2 , respectively. This shows that the above mentioned effect of the broad 
wave function leads to significant softening of the relevant scale. 

— 2 — 2 

The factorization scales Q L and Q T can be defined, in principle, in several ways. As 
a trial definition, one can look at the k 2 points, where Wi{n 2 ) reach ^Wi(0). Defined so, 
the longitudinal and transverse scales were found to be approximately equal to \{Q 2 + 
2.0 GeV 2 ) and jj(Q 2 + 2.6 GeV 2 ), respectively. The numbers g and jj are very close to 
0.15 and 0.07 0.1 obtained in 0. 

The gluon density, however, has itself significant k 2 dependence, |10j . Namely, in the 
region k 2 ~ 1 4- 10 GeV 2 and very small x g (x g < 10~ 3 ) (which corresponds, at fixed 
W = 75 GeV, to values of Q 2 < 10 GeV 2 ), Tix g , k 2 ) is a strongly rising function of k 2 . 
At larger Q 2 , the effective x g grows, and the unintegrated gluon density becomes flat. 
Finally, at large enough Q 2 (for W = 75 GeV, this corr esponds to Q 2 > 100 GeV 2 ), the 
unintegrated gluon density decreases with k 2 growth in the region k 2 ~ 1 4- 10 GeV 2 . 
Therefore, the span of effectively contributive it 2 will extend to higher values of k 2 (at 
small Q 2 ) or reduce to smaller values of k 2 (at high Q 2 ). 

In order to take this into account, it is more useful to define the hard scales via the 
following implicit relations: 

1 f°° dH 



Wi(0) 



f°° cLk 2 9 

/ — Hx g ,H)W i (H 2 ) = G(x g ,Qi); i = L,T (9) 
Jo K 
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Figure 3: The normalized weight functions Wl(k 2 ) /Wl(0) and Wt(k 2 ) /Wt(0) calculated at 
Q 2 = 100 GeV 2 . The k 2 values where Wl and Wt reach 1/2 are noticeably softer than 
(Q 2 + m 2 )/4. 
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Figure 4: The Q 2 -> % and Q 2 -> mapping in the p production. The heavy meson 
analysis expectation (Q 2 + m 2 )/^ is also shown. 

Fig. 2] shows the values of Q L (solid line) and Q T (dashed line), defined according to 
0, as functions of Q 2 . These values start from 0.63 GeV 2 and 0.4 GeV 2 , respectively, 
in the photoproduction limit, and slowly grow with Q 2 rise. At Q 2 = 27 GeV 2 (the 
highest Q 2 data point from HI data on p production), these values are only around 
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4.5 GeV 2 and 3 GeV 2 , respectively. This confirms the conclusion of |$. that even at 
largest Q 2 where data are available, we still deal with semiperturbative situation. It is 
interesting to note that a better fit to these curves is given by a non-linear, rather than 
linear approximation: 

Q 2 L « 1.5- Q 2 T ^ 0.45- (Q 2 + 1.5) - 7 , (10) 

where all quantities are expressed in GeV 2 . 

The same Figure shows also, by dotted line, the expectation (Q 2 + m 2 )/4 inspired by 
the heavy meson analysis. This expectation starts from 0.15 GeV 2 , which is noticeably 
smaller than Q^(0) and Qt(0), and rises with Q 2 significantly faster than Q L and Q T . 

3.3 The Q 2 dependence of the weight functions 

In the heavy meson approximation (JSJ), the scale Q 2 = (Q 2 + m 2 /) /A defines not only 
the hard scale of the gluon density, but also the absolute value of the cross section. It is 
interesting to check whether such a universality holds for the p production. 

This can be done by studying the Q 2 dependence of the Wl(0) and Wt(0). Motivated 
by ©, we introduce the "effective twists" e\ and e\ via the following relations: 

W L (0) = 0.153 Q • ; W T (0) = 0.153 my ■ . (11) 

\ e L> \ £ T) 

Here, factor 0.153 GeV = \/ ^ p ^o e — mp is the numerical value of all the residual 

V 3a em 

factors present in the definition of Wi. This definition of the "effective twists" e 2 and Ej> 
is such that in the heavy meson limit one would recover the familiar Q = 0.25-(Q 2 +m 2 / ). 
These "effective twists" were found to be approximately equal to 

e| ps 0.23 • (Q 2 + 1.1) ; 4 ss 0.16 • (Q 2 + 1.3) . (12) 

Again, all quantities are expressed here in GeV 2 . It is interesting to note that here, 
as well as in ()10|) . the virtuality Q 2 sums with M 2 ~ 1 -j- 1.5 GeV 2 instead of widely 
assumed m p = 0.6 GeV 2 . This is not surprising, since the mass of the p meson does not 
have much relevance to the interaction of the qq dipole with the proton. 

4 The Q 2 behavior of the p production cross sec- 
tion 

Let us turn to an issue that is closely related to the above discussion of the hard scales 
in the p production, namely, to the Q 2 behavior of the exclusive p meson production 
cross section. 

As mentioned in the Introduction, the early data on p mesons were successfully 
parametrized (in the moderate and high-Q 2 region) by a simple law 

e{l*P -» PP) oc 2 1 2 

with n = 2.32 ± 0.10 (ZEUS, 14 ) or n = 2.24 ± 0.09 (HI, However, subsequent 

data made it clear that such power fits have very limited applicability. One thus can ask 
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what would be the best way to parametrize this behavior. Note that this is not merely 
a question of how to fit the data, but rather a question of what are the main sources of 
Q 2 behavior of the cross section. 

Let us take the heavy meson limit result © as a starting point. It states that the 
cross section should quickly fall with Q 2 due to powers of Q 2 + m v , the fall being just 
slightly tamed by the rise of the gluon density. In order to check whether or not these 
are the only sources of the Q 2 behavior of the p meson cross section, let us consider a 
rescaled cross section 

{Q 2 ? • b(Q 2 ) 



Z(Q 2 ) = a( 1 *p^ pp) 



G (x g ,Q 2 ) -a s (Q 2 ) 



(13) 



(Q 2 + mi) /A. Note again that the gluon density is taken here 



Here, as before, Q 2 

at constant energy (W = 75 GeV), so that x g = 0.41 (Q 2 + m 2 )/W 2 also depends on 
Q 2 . From the heavy meson limit expression (jEJ) one might hope that £(Q 2 ) is roughly 
independent of Q 2 , since all the sources of Q 2 dependence are removed. 



rescaled cross sections 
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Figure 5: The rescaled cross section S(Q 2 ) defined according to < f?3j) . 

Fig.Elshows the ^-factorization predictions for T,(Q 2 ) as function of Q 2 . One observes 
that this rescaled cross section has quite a flat shape, from very low to very high Q 2 , 
regardless of the exact wave function used. 
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In principle, all the observables used in are accessible in experiment. The only 
delicate issue would be the choice of the gluon density, especially at low Q 2 , since the 
DGLAP fits to conventional gluon density are avaibale only for Q 2 > 1 GeV 2 and can 
differ significantly from the ^-factorization resuls, see ^Uj. It would be still interesting 
to see, whether the experimental data possess such a scaling. 

In the light of the results of the previous section, it is somewhat surprising that ()13[). 
inspired by the heavy meson limit, is nevertheless rather flat. If, however, we take a look 
at (jl and ut separately, and define 
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(14) 



we will see that this remarkable Q 2 -independence holds for S^(Q 2 ), but not for Y<t(Q 2 ), 
Fig- El left panel. This difference between Y<l(Q 2 ) and £*r(Q 2 ) can be regarded as a 
direct proof that the hard scales relevant for longitudinal and transverse cross sections 
are different. 
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Figure 6: (7e/£ 'panel) The rescaled cross sections Y,l(Q 2 ) and Ey(Q 2 ) defined according to 
P^p . T/ie Q 2 behavior of S^(Q 2 ) and E^(Q 2 ) are i>en/ different, (right panel) The rescaled 
cross sections T> L (Q 2 ) and T, T (Q 2 ) defined according to J73)). Tae two quantities are close to 
each other and both have flat Q 2 behavior. 



If we take into account the results of the previous section and introduce 
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then these rescaled quantities have flat Q 2 shape both for longitudinal and transverse 
cross sections, FigJHl right panel. Moreover, Y,l{Q 2 ) and T,t(Q 2 ) are close to each other, 
which should be expected by the construction of the "effective twists" and factorization 
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scales. These curves are not exact constants due to non-zero average value of A 2 and 
due to the presence of helicity-flip amplitudes. 

Again, it would be very interesting to see a similar analysis of the experimental data 
in terms of Xj(Q 2 ). In particular, it would be interesting to see whether such an analysis 
can help resolve the long-standing <7T-problem of the /cj-factorization predictions j!2l ll.-jj . 

5 Discussion and conclusions 

The most of the analysis of the prevous sections was done with the oscillator wave 
function. The same analysis with other wave function Ansatze will yield slightly differ- 
ent numbers in and (|12j) . Nevertheless, the general picture will remain the same. 
Namely, several observations are stable against the variations of the exact shape of the 
wave function: 

• at smaller Q 2 (Q 2 <3 -j- 5 GeV 2 ) the DGLAP factorization scale is larger than 
(Q 2 + m 2 )/4; 

• at large enough Q 2 (Q 2 > 10 GeV 2 for the transverse amplitude and Q 2 > 20 GeV 2 
for the longitudinal amplitude), the factorization scale is significantly smaller than 
(Q 2 + m 2 )/4. This should be taken as a word of caution against an unwarranted 
application of the DGLAP approach to the problem of p meson production even at 
high Q 2 ; 

• the overall Q 2 dependence of the pQCD factorization scale is significantly flatter 
than (Q 2 + m 2 )/4. This is mostly due to the specific way the /5 2 -behavior of the 
unintegrated gluon density changes, as the Q 2 increases (at fixed W). 

• the presence of m p in the often used scale (Q 2 + m 2 )/4 is misleading, since the 
p meson mass has little relevance to the color dipole interaction with the target 
proton. Instead, Q 2 appears in combinations of the form of Q 2 + M 2 with M 2 ~ 
1 ^ 1.5 GeV 2 . 

The numerically accurate understanding of the relevant hard scales in the p produc- 
tion allowed us to study in detail the Q 2 dependence of the p production cross section. 
With the help of the rescaled cross sections ()13j) - (jl5|) . Figs. |S]and|BJ we showed once 
more that the production of transverse and longitudinal vector mesons is governed by 
distinct hard scales. When we took into account the difference of the scales just found, we 
observed a very close similarity between the rescaled cross sections Tjl{Q 2 ) and T,t(Q 2 )- 
It would be very interesing to see the results of a similar analysis of the experimental 
data. This analysis might help resolve the long-standing problem of too small <jt at high 

Q 2 . 

Let us also comment on a contribution to the recent puzzle of energy independence 
of r p = a(j*p — » pp) / <Ttot{l* p) ratio. As was discussed in this paper, at higher Q 2 , the 
true factorization scale is smaller than (Q 2 + m 2 )/4. Inverting the argument, one can 
state that <?tot{Q 2 t = (Q 2 + m p)/4) is expected to be harder than a p {Q 2 ). This should 
enhance the expected value of the Pomeron intercept in otot and reduce the prediction of 
the energy behavior of r p . Unfortunately, our numerical analysis showed that this effect 
is marginal and does not lead to resolution of the problem. 
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